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A new technique named Generalized Borel Transform (GBT) is applied to the generating 
functional of the <E> theory in zero dimensions with degenerate minima. The analytical solution 
of this function, obtained in the non perturbative regime, is compared with those estimations 
predicted by large order perturbation theory It was established that the GBT is a very efficient 
technique to capture these contributions. On the other hand, renormalons associated to the 
resummation of those perturbative series were not found to be the genuine source of the non 
perturbative contributions of this model. 
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Substantial improvement of experimental data has recently demanded theorists to give 
more accurate estimates since the finite order perturbative predictions have uncertainties 
comparable with experimental errors. This has renewed the interest in computing non 
perturbative contributions from the resummation of series. In Quantum Field theory non 
summable Borel, the problem of obtaining these contributions has been largely studied in the 
literature. 1 Among the different methods proposed to tackle this problem, the Conventional 
Borel Transform (CBT) of the weak coupling expansion 2 and its modified version (MBT) 3 ' 4 
have become popular tools for the reconstruction of the physical quantities predicted by 
these theories. Basically, these techniques require the exact knowledge of the large order 
behavior of the perturbative series in the coupling parameter. In spite of the lack of these 
series, which are divergent, 5 one can make use of their coefficients to define other expansion 
in powers of a new variable named Borel and in such a way that it is convergent. Then, the 
resummation of the initial asymptotic series is formally recovered by integrating the con- 
vergent series on the Borel variable. In this approach, the non perturbative contributions 
are estimated from the ambiguity appearing in the integration, where singularities (poles 
and/or branches) called IR renormalons 6 localized on the positive real semiaxis of the Borel 
complex plane invalidate these proposals. On the other hand, an alternative development, 
called the Generalized Borel Transform (GBT), 7-9 was recently presented which has proved 
to be potentially useful in overcoming those difficulties. It avoids the use of series by intro- 
ducing auxiliary parameters having no physical meaning which are introduced for the sole 
purpose of making the calculation analytically tractable. These parameters together with 
the analytic properties of the GBT are the fundamental tools needed to obtain analytical 
solutions of parametric integrals like Laplace-Fourier-Mellin transforms (LFMT) for all the 
range of its parameter. Therefore, it is extremely useful to study non-perturbative regimes. 
In addition, it involves simple mathematics which requires basic calculus such as the evalu- 
ation of derivative, indefinite integrals and limits. In fact, this technique has already been 
successfully applied to different areas of Physics such as Quantum Field Theory, 7 Polymer 
Physics 8 and Quantum Mechanics. 9 

In connection with the techniques previously mentioned, it is certainly instructive to 
study exactly solvable models so that the efficiency and precision of the different proposals 
can be quantitatively checked. This analysis may help in defining criteria for selecting the 
appropriate approach to be used in more realistic models. In this sense, it is largely accepted 
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that the analysis of zero-dimensional models can provide some insight on the behavior of 
Green's functions in Field theories and Statistical systems. 10 In particular, the understand- 
ing of some aspects of zero-dimensional Field Theories could leave trails about the non 
perturbative behavior of these theories in higher dimensions. For instance, classical actions 
with degenerate isolated minima play a fundamental role in studying processes governing 
the decay of metastable atomic and nuclear states as well as the transition of overheated or 
undercooled thermodynamic phases to a stable equilibrium phase. 

One of the most popular model of this kind is the well-known zero-dimensional <3> 4 theory 
with a degenerate minima 1 due to its purely non perturbative feature. The analysis of this 
model is often based on the study of its generating function, starting point for computing 
moments and Green's functions. Moreover, at present the exact solution is known from the 
numerical computation of its integral representation 

/+oo 
dxexp[-x 2 (l-gxf /2] , g > 0. (1) 
-oo 

Therefore, the main purpose of this paper consists in discussing the use and to remark 
the advantages of the different approximate computing tools aforementioned to obtain an 
analytical approximate solution of the generating function JH) in the non perturbative regime. 

Thus, the analysis of this model is initially done on the context of Perturbative Field 
Theory in which non perturbative contributions are estimated from the solution of the 
generating functional JQ) in powers of g. It is easily obtained by rewriting the integrand in 
Eq. dH) as a product of two exponential functions as follows 

/+oo 
dx exp [-x 2 /2] exp [gx 3 - g 2 x 4 /2] . 
-oo 

Then, the expansion of the second exponential in powers of its argument leads to the 
following formal solution of this model 

Z P (g) = v^F + V2± F < 2 " + 5/2) (8 g y +1 , 9 ^0+. (2) 

n=0 ^ ' 

Even though this solution is exact, it is plagued of difficulties and consequently it does 
not represent a suitable analytical solution. Indeed, the series is asymptotic and factorial 
divergent. Basically, it means that the contributions of the terms of the series decrease with 
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increasing n for n < n m « n and increase for n > n m i n . The value of n m i n = N Q (g) corresponds 
to the term with the minimal contribution to the series. Thus, the best approximation 
provided by expression Q is obtained by truncating the series at n min . Given that it 
depends on the value of the parameter g, then the accuracy of the result will also depend 
on it. For example, small values of g provide a large N Q (g) and consequently there are 
sufficient amounts of terms contained in this series so that a good approximation of the 
exact solution is obtained. In contrast to this, values of the parameter g going away from 
the perturbative regime reduce the value of N Q (g) and the series involves just a few terms. 
Thus, it is significantly off from the exact one. In addition to this trouble, their coefficients 
do not alternate in sign avoiding the cancelation of the consecutive terms in the series, which 
makes it worse the divergent behavior. 

In some cases, the aforementioned difficulties can be minimized, or hopefully eliminated, 
by applying standard approaches specifically designed to deal with this class of series. In this 
model, it is clear that the use of the exact expression of their coefficients would simply recover 
the starting expression for the generating function given by the integral representation JTJ) 
whose analytical solution is unknown. However, an approximate analytical solution for this 
function can be obtained from Eq. J2|) by using the following asymptotic behavior for their 
coefficients 

in which case the resulting approximate generating function becomes a factorial divergent 
series 

^(^) = v / 2^+^^r(n + l)(32^ 2 ) n+1 , g^0 + . (4) 

which can be resumed by applying the Conventional Borel Transform as follows. The 
Borel sum of the generating function is obtained from expression (jlj) by using the integral 
representation of the Gamma function. 11 Afterwards, the sum and the integral operations 
presented there are formally interchanged to get 

r+oo 

Z Borel (g) = V2^+ dt exp [-t] B (t) , (5) 
Jo 
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where B (£) is the Borel Transform defined by 

Since the last expression is valid only for "52tg 2 < 1, the integral in Eq.(|s3) can be com- 
puted whenever an analytical continuation of its integrand to the Borel complex plane is 
available. In doing so, it appears a singularity at t = l/32g 2 . As a consequence, the resum- 
mation becomes invalid (non summable Borel). Anyway, a meaning to the ill-defined integral 
(jSJ is often given by adopting a regularization prescription. Nevertheless, this generates an 
ambiguity in its evaluation since the result will depend on the particular prescription taken. 
If the Cauchy Principal value is used, the pole is avoided by deforming the contour of inte- 
gration along two semicircles centered at the pole above and below the real axis respectively. 
Hence, 6 ' 12 



. LU!rd (g) = V^ + ^PV / +0 ° dt ^ V \J\ ± zv^exp [-l/32<? 2 ] (7) 
V?r J 1 - 32tg 2 



where the last term comes from the residue evaluated at the pole. 

From this definition, the estimation of the non perturbative contributions is immediate. 
The generating function defined initially by equation JTJ is analytic in the whole complex 
plane of the coupling parameter g. Then, the non perturbative contributions which had been 
excluded in the computation, must cancel those non analytical contributions present in the 
expression in order to preserve this property. In this sense, the non analytical structure 
of the Borel transform is assumed to contain information about the origin of the main non 
perturbative contributions of the model, leaving aside the ones unrelated to renormalons. 

Therefore, within this technique one ends with an approximate solution for the generating 
function in the form 3 



Z Borel (g) = V 27T H — exp 
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^(32^) +CeXp[_1/32 ^ (8) 



32<? 2 . 

where Ei (z) is the exponential integral function 11 and C is a constant to be adjusted from 
"experiment". In this model it means from the numerical exact solution (fT|. 

The numerical evaluation of this result is presented in Fig. 1. This plot clearly shows 
that the prediction coming from the large order perturbative theory is not good for the non 



5 



perturbative regime. Therefore, the non analyticity of the Borel sum (IR renormalons) is 
not the true source of the main non perturbative contributions of this model. 

In more realistic theories the exact expression of the coefficients of the series is otherwise 
rarely known. 1 However, one can analyze other asymptotic expansion of the coefficients in 
this model which converge more rapidly to the exact one. Concretely, one of them reads 



Z P (g) = V2tt + — ^ , g -+ . (9) 

V 71 n=0 (n + 1) 

By applying the previously explained Borel approach to this series, the Borel Transform 
adopts the form 

r i/2 - (32^ 2 ) n , . 

^(^ + 1) 1/2 



This series converges to the Polylog function 11 L 1 / 2 (32tg 2 ) in the interval [0, l/32g 2 ] . 
Explicitly, the expression (jHJ) becomes 



J Borel 



i r+oo 

(g) = v / 2^+ —j= \ dtexp [-t] t~ 1/2 L 1/2 (32tg 2 ) (11) 
v n Jo 

where right now the integrand in Eq.JHJ) must be analytically continuated to |32tg 2 | > 1. As 
a consequence, it acquires a cut on the positive real semiaxis in the interval [l/32g 2 , oo] due 
to the well-known Polylog function analytical properties. 11 As expected, the resummation 
becomes invalid again. In this case, the non analyticity is avoided by deforming the contour 
of integration along two parallel lines slightly above and below the real axis. Hence, it 
reads 12 



Y r-roo 

Z B orei {g) = V2^+^ dt exp [-t] t~ 1/2 Re {L 1/2 (32tg 2 ) } (12) 
V n Jo 



with 



+oo 

dtexp [-t] r 1/2 Im {Li/2 (32tg 2 )} 

l/32 ff 2+ 



Im{L 1/2 (32tg 2 )} 



7T 



In (32tg 2 
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being the non analytical part giving an estimation of the non perturbative contributions. 
Therefore, the suitable approximate solution for the generating function has the form 3 



r+oo 

Z Borel (g) = — dt exp [-t] t- l ' 2 Re {L 1/2 (32tg 2 ) } (13) 
V 71 Jo 

+c r 

il/32 9 2 V ln ( 32 V) 

Unfortunately, the analytical solution of these integrals is unknown. The numerical eval- 
uation of this result (Fig. 1) shows that the approximation (THfl) improves the prediction of 
non perturbative contributions when it is compared with the first approach (jHJ). In spite 
of this improvement, the expression (fl~3|) still does not provide a good estimation of the 
generating function (QJ in the non perturbative regime. 

Setting aside the difficulties present in the Borel resummation technique, the computation 
of the generating function by using the GBT is now presented. In doing so, it is necessary to 
rewrite the expression (0) in terms of a Laplace transform. It is achieved as follows. Firstly, 
the variable of integration x is replaced by u/g and a new coupling parameter a = 1/ [2g 2 ] 
is defined. Secondly, the integration on the whole real axis is performed on each semiaxis, 
positive and negative, separately. In addition, the last integral is rewritten on the positive 
semiaxis by changing variables u = —v. As a consequence, the expression (JTJ) becomes 



r r+oo r+oo \ 

Z (a) = \pla < / duexp \—au 2 (1 + w) 2 ] + / dvexp [— av 2 (1 — v ) 2 ] > . 



(14) 



Then, the second integral in Eq. lfTHl is divided into plus J^°° so that the change of 
variables w = 1 — v transforms J^°° into the first integral in Eq. |T4| . Finally, this is written 
as a Laplace transform by changing variables y = u 2 (1 + u) 2 to get 



Z{a) = V2a~{ dy . y> =+ / dv exp \-av 2 (1 - v) 2 } } . (15) 
[Jo y/y^l + A^y Jo J 

Observe that this new expression for the generating function presents important advan- 
tages. A simple analysis of the behavior of both integrals in Eq. (|T5l) allows one to identify 
what is the source of the non perturbative contributions. In fact, they show that, in the 
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limit g — > oo + (a — > + ) , the main non perturbative contribution comes from the Laplace 
transform. Furthermore, it can be analytically computed by applying the GBT technique 
which is described below. 

Basically this method consists of introducing two auxiliary functions, S (a, n) and 
B\(s, n)(the Generalized Borel Transform) which depend on auxiliary parameters called 
n and A. These parameters have no physical meaning and are introduced for the sole pur- 
pose of helping in the computation of an explicit mathematical expression for the following 
Laplace transform 

r+oo 

S (a) = / duH (u) exp (—au) , a > 0. (16) 
Jo 

In doing so, S (a, n) is defined in terms of S (a) by the formula 

an poo 

S{a,n) = {-) n —S(a)= u n H (u) exp (-au) du, n>0. (17) 
oa n Jq 

which can be inverted to get 

/„ n— 1 

da-- - j da S (a,n) + ^^c p (n) aP . (18) 
» ' p=° 

Note that the finite sum comes from the indefinite integrations and it can be omitted 
whenever the Laplace transform (|T6|) fulfills the following asymptotic condition 

lim £ (a) = 0. (19) 

a— >oo 

In addition, the approximate computation of S(a,n), and subsequently of S (a) from 
Eq. (fT8|) . is doable from the following definition of the Generalized Borel Transform of S (a, n) 



B x (s, n) = - 



exp [s/rj (a)] 



-i -As 



A77 (a) 



S (a, n) drj (a) 
lv(a)} 2 da 



da, Re (s) < 



(20) 



where A is any real and positive number, whereas rj is defined so that \jr\ = A (exp (a/A) — 1). 
Then, the analytical properties of B\(s,n) allow one to invert unambiguously Eq. ([20"|) to 
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get 



OO POO 



S(a,n) = 2X (1 — exp (— a/A)) / / exp [G (w, t, a, A, n)] dwdt. (21) 



— oo J — oo 



being G (w, t, a, A, n) an involved function depending explicitly on if [8, 9]. 

Given that the aforementioned expressions (|17II21|) are valid for any number of n > 0, the 
dominant contribution to the double integral is computed for n 3> 1 by using the steepest 
descent method 13 in the variables t and w. On the other hand, observe that each value 
of the parameter A in Eq. (|2T)l) defines a particular Borel transform. However, the resulting 
expression for S(a,n) given by Eq. (|2T| does not depend on A explicitly. Consequently, one 
can choose the value of this parameter in such a way that it allows one to solve Eq. (|2Tf . 
Hence, the approximate expression of S (a, n) obtained from the saddle point is computed 
in the limit A — > oo to finally obtain 

Sa P tox (a, n) = V2ire~ 1/2 )/JJ?tl + 1 [u ] n+1 H (u ) exp [-/ (u )\ , (22) 

V D ( u o) 

where 

/ (u ) = l + n + Uo (23) 

D K) = -u ^p^- [1/2 + / (u )} + f (u ) [1 + f (u )} , (24) 
du a 

and u Q is the solution of the following implicit equation 



u 2 a 2 = f (u ) [f (u ) + 1] . (25) 

In summary, the application of this technique consists simply in solving the implicit 
equation (|25|) for n 3> 1 to obtain the expression of u Q and replace it in Eq. ([22|) . Then the 
n— indefinite integrals appearing into Eq. (fT8|) can be solved to get an approximate solution 
of S (a) in the limit n — > oo. Observe that the expressions (jl8H25jl appearing in this approach 
avoid perturbative expansions and all of them are unambiguously defined. Moreover, the 
parameters n and A initially introduced to make the computation mathematically tractable, 
are lastly taken away from the approximate solution for S (a) . 
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Thus, the approximate expression for the Laplace transform appearing in Eq. (|15|) is 

into the expressions 



obtained by replacing the function H (u) = 1/ y/u\J\ + k\Ju 
provided by the GBT. Then, the expression of / (u ) given by expression (f23|) is replaced 
into Eq. ([2~5T) and n is assumed large to get u Q ~ (n + 3/4) /a. Then / (u ) — n + 3/4 , 
D (u ) ~ (n + 1/4) (n + 5/4) and S (a, n) can be approximated as 



2irJn + 5/4 [n + 3/4] n+1/4 exp [- (n + 3/4)] 



2y/(n + 1/4) (n + 5/4)a n+1 /\ 1 + y/aj 4^ + 3/4 



^Apr-ox (O; H.) 



By recognizing the asymptotic expansion of the Gamma function 11 r (n + 1/4) in Eq.(|2 
the expression of Sa P tox can be simplified to get 



(26) 



S 



Aprox 



(a, n) 



T(n + 1/4) 



2a n+1 / 4 Wl + y/aj Ayjn + 3/4 



(27) 



Then, this expression is replaced into Eq. (|18|) . and the coming expansion 



1 + vV 4^n + 3/4 p 



-1/2 



(0) 



7 p/a 



^ (-4f („ + 3/4f / 2 



(28) 



is utilized, being (x) the Laguerre polynomial. This leads to calculate the non pertur- 
bative contribution of S (a) as follows 



Snp (a) ~ 



lim 



n — > oo 



^ I> + 1/4)L; 1/2 (0) 
^ 2 (-4f („ + 3/4f / 2 



da- ■ ■ / da 



1 



n+l/4-p/2 



(7 



(29) 



where 



da - ■ ■ / da 



1 



a 



n+l/4-p/2 



-) n r(l/4-p/2)/[r(n + l/4-p/2) 



1/4-P/2] 



Thus, the expression (|29|) computed in the limit n — > oo turns into 
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^ L p 1/2 (0)r(l/4-p/2) 

-W (a) - 2^ 2F4)W ' 

p=0 v 7 

whose series converges to the following final expression 

s -w-HMI) + MI0M-|)> < 30 > 

where (2) is the Bessel function of second class. 11 

Going back to expression (fl~5|) . the solution for the generating function is completed by 
solving the second integral appearing there. It is easily done by using the expansion of the 
exponential term and exchanging sum and integral operations. Then, the following solution, 
valid for any value of the parameter a, is obtained 



„1 / 

/ dxexp \— ax 2 (1 — x) 1 = > 



-a) k [T (2 + 2k)f 



- (2k + iyT(k + l)T(2 + 2k) 



2^2 









"3 5" 


( 


Ik 




4' i 



a 

~16 
(31) 



where 2-^2 is the well-known Generalized Hypergeometric function. 11 

Finally, the expressions (f30|) and (|3T| evaluated in terms of the original coupling pa- 
rameter g provide the following analytical solution of the generating function in the non 
perturbative regime 



(9) 



'1/4 



64# 5 



1-1/4 



exp 



64# : 



(32) 









"3 5" 




( 






.4' 4 


' 32^)] 



+2F2 

The test of the accuracy provided by the GBT result is presented in Fig. 1. The relative 
error with respect to the exact solution (JTJ (solid line) is less than 5% for values of g 
larger than 100. It shows that the GBT is a very efficient tool to capture non perturbative 
contributions of this model. 

In summary, expressions (JB|), (|T3|) and (|32|) for the generating functional of the zero- 
dimensional $ 4 theory with degenerate minima in the non perturbative regime of the cou- 
pling parameter g were obtained by application of the Borel resummation and GBT tech- 
niques. Models like this are mainly utilized for studying the physics of barrer penetration 
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and tunelling effects in systems which classical minima are connected so that the symmetry 
between them is not spontaneously broken. 

It was found that the genuine origin of the main non perturbative contributions of this 
model is not the Borel ambiguity. This non analyticity of the Borel sum is due to the 
non alternating character of the coefficients in the expansions in powers of g. In fact, this 
behavior is a direct consequence of the degeneracy of the ground state and the presence of 
inevitable contributions, which cannot be captured from Perturbation Theory. It implies 
that for reconstructing the generating function one needs additional information which is 
absent in Perturbation Field Theory. Indeed, this approach has to introduce a parameter 
to be afterward fitted in order to quantify non perturbative contributions. 

On the other hand, the GBT was established to be a very efficient tool to capture non 
perturbative contributions of this model. The overall goal of this technique comes from the 
fact that it involves expressions unambiguously defined which preserve the complete infor- 
mation about the coupling parameter. Indeed, the computation of the analytical solution 
is doable due to the help of auxiliary parameters which do not appear in the final results. 
The main role of the real parameter of the model is simply to control the localization of the 
saddle point u Q and consequently the behavior of the resulting final solution. 

The author thanks M. Pettitt and the Institute for Molecular Design, Houston for the 
generous and warm hospitality extended to him. He also thankful H. Fanchiotti, Luis Epele 
and C.A. Garcia Canal for their advice and helpful discussions. 
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FIG. 1: Generating function Z (g) as a function of the coupling parameter g. The solid line 
corresponds to the numerical evaluation of the exact solution (QJ and the dashed line to the GBT 
solution (JS21- Dotted-dashed and dotted lines correspond to Borel resummation solutions (fT3")l and 
|(EJ respectively. The value of C minimizing the deviation of these approaches with respect to the 
exact solution Q (solid line) was found to be equal to —47.313 and 3.572 respectively. They were 
computed by the least square method in the range of values of g given by the interval [100, 1000]. 
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